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Abstract 

The cross sections of e + e~ to S— wave and P— wave quarkonia with C-parity even associated with 
a photon are analyzed in the framework of non-relativistic Quantum Chromodynamics(NRQCD) 
factorization formulism. The short-distance coefficients are analytically determined up to the next- 
to-leading order(NLO) QCD radiative corrections to S— wave and P— wave quarkonium production 
and NLO relativistic correction to r\ c production. We also invoke the analytical expressions to 
estimate the cross sections. Our numerical results indicate that both the QCD and the relativistic 
corrections are considerable. 

PACS numbers: 12.38.-t, 12.39. St, 13.66.Bc, 14.40.Gx 
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I. INTRODUCTION 

Non-relativistic quantum chromodynamics(NRQCD) factorization formalism [l| is an use- 
ful tool for analyzing the inclusive production of the heavy quarkonium. According to it, the 
cross section is expressed as a sum of products of short-distance coefficients and NRQCD 
matix elements. The short- distance coefficients can be calculated as perturbation series in 
coupling constant a s at the scale of the heavy quark mass. The matrix elements scale in 
a definite way with the typical relative velocity v of heavy quark in the quarkonium state. 
Thus the production cross sections can be expressed as double expansions both in a s and 
in v to any desired order. 

Recently, the cross section of the exclusive process e + e~ —> H + 7 at the center-of- 
momentum(CM) energy y/s = 10.58GeV was computed by Chung, Lee, and Yu^, where 
H represents a S— wave or P— wave heavy quarkonium state with charge- conjugation parity 
C = +1, to the leading order (LO) both in a s and in v in the NRQCD factorization formalism. 
Their predicted cross sections are about 82fb, 49fb, lfb, 14fb and 5fb for r) c , rj c (2S), Xco, Xa 
and Xc2 production, respectively. Some early calculations of similar processes like Z° — > H+~f 
were presented in [3| in the color-singlet model. 

These relatively large cross sections make the processes be measurable at the B factories. 
To have more precise predictions, it is necessary to carry out the calculations to higher 
order both in a s and in v. One may expect that those correction terms are considerable. 
Especially, for charmonium production, both a s and v 2 are not small. In this paper, in the 
framework of NRQCD factorization formalism, we compute the short- distance coefficients 
of the leading NRQCD matrix element to the NLO QCD corrections both for S— wave and 
P— wave heavy quarkonium production. For the r) c , we also calculate the tree level short- 
distance coefficient of the v 2 relativistic correction term. We will not perform the calculation 
for the NLO relativistic corrections to t] c (2S) and P— wave quarkonium production since the 
NLO NRQCD matrix elements for these states are still not determined preciously by fitting 
data. The NLO QCD correction to the e + e~ — > r] c + 7 was calculated first in J4]. However, 
we found that their original expressions given in Eqs. (23), (25) in their paper were incorrect. 
Probably there were some typos in their equations. 

The NRQCD factorization formulism separates the short- distance effects happened at 
energy scale m or higher, where m is the heavy quark mass, and the long-distance ef- 
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fects happened at long-distance 1/mv or larger which describes the formation of the heavy 
quarkonium. The short-distance effects are described by the short-distance coefficients while 
the long-distance ones are described by the NRQCD long-distance matrix elements. In the 
process e + e~ — > H + 7, the emission of the photon can be either short- distance effects or 
long-distance ones depending on the emitted photon being hard or soft comparing to the 
quark mass m. Thus when the NRQCD factorization approach is used to analyze the process 
of e + e~ H + 7 at the B factories, the factorization formulism for the cross sections takes 
different forms in the cases of H being cc or bb state. In the former case, the emitted photon 
is so hard that the short-distance process is e + e~ — > cc + 7. However, in the later one, the 
energy of the emitted photon is much less than the b quark mass so that the emission of the 
photon is a long-distance effect and should be described by the long-distance matrix as an 
electric El or a magnetic Ml transition. 

The rest of this article is organized as follows. In Section II, we apply the NRQCD 
factorization formulism to analyze e + e~ — >■ H + 7 at the B factories where we distinguish 
two different cases of cc or bb, corresponding to the emitted photon being hard or soft. In 
Section III, we determine the short-distance coefficient for the cross section of the process 
e + e~ — > f] c + 7 up to the NLO radiative QCD correction for the LO matrix element and the 
LO short-distance coefficient for the relativistic correction term. We also discuss various lim- 
itations of those short-distance coefficients. In Section IV, we determine the short- distance 
coefficients for the cross sections of the processes of e + e~ — > Xcj + 7 up to the NLO QCD 
radiative corrections and discuss various limitations. Section V contributes to discussions 
of the factorization form for the processes of e + e~ — > ^(xbj) + 7 with the emitted photon 
being soft. In Section VI, we apply the obtained formulae to carry out numerical estimations 
for the cross sections of those processes. We also summarize our results in this section. We 
present the analytical NLO short- distance coefficients of the leading order NRQCD matrix 
elements in v for the cross sections of e + e~ — > r] c + 7 and e + e~ — > Xcj + 7 in Appendix A 
and Appendix B, respectively. 

II. NRQCD FACTORIZATION FORMULISM FOR e+e H + 7 

In the process of e + e~ — > H + 7, where H represents a S— wave or P— wave heavy 
quarkonium state with C-parity even, an almost on-shell and near threshold QQ pair is 
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created at the short-distance scale of 1/m or smaller, then constitutes a heavy quarkonium 
H that happened at the distance scale of 1/mv. The production cross section can be analyzed 
in the NRQCD factorization formulism |l| as double expansions both in velocity v and in 
coupling a s . The factorization takes two different forms in two different cases corresponding 
to the emitted photon being hard and soft. 

For charmonium production at the B factories via process e + e~ — > charmonium + 7, 
the momentum of the emitted photon satisfies \k\ > m. The photon is so hard that the 
short-distance effects arise from the process e + e~ — > cc + 7. We assign K. n with n in the 
collection of 1 S , 3 Po, 3 Pi and 3 P 2 as follows [lj]: 

/C 3Pl = -j={-\%> x °")> 

/Ca P2 = - Z -D^\ (1) 



The NRQCD factorization formulism for the cross section of e + e — > H + 'jhard (H = r], 
r] c (2S), Xcj) takes the following uniform form: 



CI 



+ ^URe(0|x t /C^|^)(^|^X(-^D) 2 X |0), (2) 

where (0|x f /C n ^|i/) (H\^)C n x\0) a,ndRe(0\ X j IC n ilj\H}(H\^]C n (-^B) 2 x\0} are NRQCD ma- 
trix elements and Fi(n) and Gi(n) are the short- distance coefficients corresponding to the 
matrix elements. The leading contributions to the matrix elements arise from these states, 
which possess the same quantum numbers as n. In ([2]), the NRQCD matrix element in the 
second term is suppressed by v 2 compared to that in the first term. The second term is usu- 
ally called as the relativistic correction term. The short-distance coefficients can be expanded 
as power series of a s at energy scale of m or higher. They can be determined by matching 
the cross section of a free quark pair production process e + (p e i) + e~(p e 2) — > c{pi)c{j>2) + l{k) 
with the color-singlet on-shell cc pair near the threshold. Given the total and the relative 
momenta of the cc pair are P and 2q, respectively, we then have 

Pi = -P + q, P2 = ^P~q, pl=pl = m 2 , P-q = 0. 



Then the S-matrix element for this process is given by 



iM 



(3) 



where e is the electromagnetic coupling, e is the polarization four-vector of the photon and 
s is the square of the center-of-mass(CM) frame energy. The leptonic current L M in ([3]) is 
defined by 



Lf, = v(p e2 )7^(Pei), 



and the amplitude A 11 ' is given by 



A» v = Tr [A^Ur 



(4) 



(5) 



where IL m is the projector operator of spin-single (m = 0) or spin-triplet (m = 1) state and 
A 1 " is the amplitude of e + + e~ — > cc + 7 with the wave-functions of the external quark 



lines removed. For the spin-singlet state, n reads 



5] 
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(f/2- 4-m) lb (f + 2E){r/2+ 4+m) 



+ m) 

and for the spin-triplet state with spin-polarization vector e, ITi readsj^] 

(f/2- i-m) f(P + 2E)(P/2+ 4+m) 



IT 



+ m) 



where E = \Jm 2 + q 2 . 



(6) 



(7) 
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FIG. 1: The tree level Feynman diagrams for e + e — > cc + 7. 



At tree level, there are two Feynman diagrams contributing to the process as illustrated 
in Fig. [TJ A^ u is written as 



A" v = e%e 2 



m 



V 



(p2 + k) 2 — m- 



u . „ Pi 
7 +7 



m 



(pi + k) 2 — m 



(8) 
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where eQ is the electric charge of the heavy quark. 

By expanding the amplitude in terms of q, the contributions arising from S— wave and 
P— wave can be picked out. The cross section in the QCD can then be evaluated by squaring 
the amplitude and integrating over the phase space. The cross section of the same process 
can also be evaluated in the NRQCD factorization formulism with the hadron in the NRQCD 
matrix elements replaced with the free quark pair being of the same quantum number as 
H. The short- distance coefficients Fi(n) and Gi(n) can then be inferred from equaling the 
cross section calculated both in the QCD and in the NRQCD factorization formulism. In 
our paper, we determine the Fi( l So) and Fi( 3 Pj) to the NLO QCD corrections in a s and 
Gi( 1 S'o) to the LO in a s . 

On the other hand, for bottomonium production at the B factories via process e + e~ — > 
bottomium +7, the momentum of the emitted photon yields |A;| <C m. The photon is so soft 
that the wave-length of the emitted photon is compatible with the size of the bottomonium. 
Hence the emission process is sensitive to the long-distance effects. Its effects should be 
absorbed into the long-distance NRQCD matrix elements. The short-distance effects arise 
from the contributions of the process e + e~ — > bb with bb in a spin-triplet and a color-singlet 
state. So, the NRQCD factorization formulism for the cross section of e + e~ — > H + j so ft 
takes the following form: 



* = ^^(0|xVV|^ + 7)(^ + 7l^Vx|0) 



A t^-S/,| U _L rJ\lW _L /vU/it/T*"-, 

+ ^!^ Re <o| x t (7 V|tf + 7><tf + Yl^(-^D)Vx|0>, (9) 



where -F 2 ( 3 Si) an d G2( 3 Si) are the short-distance coefficients. Again, those short-distance 
coefficients can be inferred from perturbative QCD. The long-distance matrix elements are 
non-perturbative in nature. 

The short-distance coefficients in this case can be determined by matching the cross 
section of the free quark pair production process e + + e~ — > bb( 3 Si), with the bb pair near 
the threshold. The cross section of this process can be evaluated both in the QCD and 
in the NRQCD factorization formulism (Q. The short-distance coefficients can then be 
determined by equaling the cross section inferred from both theories. The short- distance 
coefficients -F 2 ( 3 S'i) and G , 2( 3 5'i) in © have been determined to the NNLO0[3] and to the 



NLO 



8| in the expansion of a s in literatures. 
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III. e+e -> 7y c + 7ft, ar . d 

In this section, we apply NRQCD factorization formulism (J2J) to the process e + e _ — > 
Vc + lhard and determine the short-distance coefficients F 1 ( 1 S , ) and Gi( 1 5 , ) to the NLO and 
to the LO in a s respectively. As mentioned in previous section, they can be determined by 
matching the process e + (p e i) + er{p e 2) —> c(pi)c(P2)( l So) + n f(k) with cc in the color-singlet 
1 5'o state near the threshold. 

Since the cc is in the 1 S'o state, Lorentz and CPT invariance implies that the tensor A^ v 
can only be expressed as 

A» v = Ae^P^ , (10) 

where A is a Lorentz invariant quantity and can be inferred from perturbative QCD. Con- 
sequently, the differential and the total cross section read: 

da a\A\ 2 , AE 2 



3 a + x 2 ), (ii) 



dx 64 s 
a\A\ 2 , AE 2 



Here, we use a notation x = cos#, where 9 is the angle between the photon and the beam 
line at CM frame. 

The cross section for the same process can also be derived in the NRQCD factorization 
formulism. Thus, we can determine the short-distance coefficients by equaling both results. 



A. Short-distance coefficients in the LO in a s 

In this subsection, we determine the short-distance coefficients F 1 ( 1 S , ) and Gi( 1 S , ) to 
the LO in a s . They can be obtained by matching the cross section of e + e~ — > QQ( l So) + 7 
at tree level both in the QCD and in the NRQCD factorization formulism. To gain G 1 ( 1 S'o) 
which is the short-distance coefficient of the relativistic correction matrix element, we need 
to expand the cross section of this process to order q 2 . 

Imposing the on-shell condition of the external quarks on (jSJ), A^ v in (jSJ) is simplified as 

A» v = e 2 Q e 2 Tr 

The Lorentz structure of the projection operator Ho can be written in a general form 

n = a 75 P + b l5 (P 4-4f>) + c 7 5 , (14) 
7 



7 V fcrf 7^ fcY 

2pi ■ k 2p 2 ■ k 



IT 



(13) 



where a, b, c are Lorentz invariant and can be easily inferred from ([6]). One see that only 
the first term in (fT4"j) gives contribution to (fl"3]) . Inserting this term into ffl"3~l) . we find 



^ = -6*ae 2 e 2 ( — ^- + — ^- ) e^P a k^ . (15) 



It readily follows that 



-6iae 2 e 2 P-k -24iae 2 e 2 / q 2 \ q 4 

where we have performed a substitution 

<?<t ~> o + K • ( 17 ) 



3 V 4E 2 
Meanwhile, expand a with respect to q 2 , to obtain 
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Combining (fT6|) and (fl"8j) . we present the total cross section (fl2l) in the expansion of q : 
as: 

CT (o) = _^ (1 _ r) 

s z 

64e> 2 (l + 2r)a 3 q 2 

where r = Am 2 /s and superscript (0) denotes the contribution arising from the LO in a s . 

Next, we turn to calculate the production cross section of e + e~ — > cc^Sq) + 'jhard in the 
NRQCD factorization formulism (J2J). At tree level, the matrix elements give: 

(cc( 1 5 )|^ t x|0) = 2Ev / 6, 
(ccCS )\^(- l -B) 2 X \0) = 2EV6 q 2 . (20) 

Inserting them into ([2]) and then expanding (2) in terms of q 2 up to order q 2 , we obtain 

a (o) = 24F 1 (0) ( 1 5 )(1 + 4) + 24Gi 0) (^o)4 ■ ( 21 ) 

By equaling (1211) with (119j) . we determine the short distance coefficients Fi^Sq) and Gi^Sq) 
to the LO in a s : 



32e%n 2 a 3 



e? ) ( 1 So) = - ~-J s2 - (i-^. 
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Finally, the cross section of e + e — > 7] c + 7 to the LO in «<, and the NLO in v 2 is given as: 
r(0) = 8e q 7ra (1 -rXOlxVl^cX^l^xlO) 

i «-* _ 

(23) 



3m 2 s 



2 

4 ^2^,3 



32envr a 

Q ^ 4g2 (1 - i )Re[(0|xVk c )(r/ c |^ t (-2 D ) 2 x|0)]. 



Here the two hadron matrix elements are normalized relativistically, which is related to the 
standard non-relativistic normalized hadron matrix elements lOf] by 



\XHW) = V^mX^\Vc)BBL + ^M(0\xK~V) 2 ^\Vc)BBL]), 



1 a 



(oixK—W^Vc) = V±m(0\xH-~B) 2 i>\i lc ) BBL . 



(24) 



Inserting ( 124|) into ( 1231) . we obtain the cross section with the hadron matrix elements in 
non-relativistic normalization as 



a 



(0) 



32e 4 Q n 2 a 3 

3ms 2 
16eg7r 2 a 3 

9m 3 s 2 



l-r)(0|xVl^)(r/ c |^x|0) 



BBL 



l its 



-(5 + r)Re[(0|xVkc>(^|^(-oD) 2 x|0) 



BBL\ 



(25) 



B. The NLO QCD correction 





FIG. 2: One- loop Feynman diagrams for e + e — )• cc + 7. 
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In this subsection, we determine the short-distance coefficients Fi( 1 S'o) to the NLO in 
a s . To this end, we need to match the cross section of the process e + e~ — > cc(}So) + ^hard 
to the one-loop level but the LO in v 2 both in the QCD and in the NRQCD factorization 
formulism. 

The Feynman diagrams responsible for the one-loop QCD correction are illustrated in 
Fig. [2l We perform the renormalization in on-mass-shell(OS) scheme and take dimensional 
regularization to regulate both the ultraviolet (UV) and the infrared(IR) divergences. In OS 
scheme, external quark lines do not receive any QCD corrections, and the counterterms from 



the renorma 
are given by 



ization constants of the heavy quark wavefunction and the heavy quark mass 



Hi: 

6Z™ = -C^f^- + A_3 7s + 31n^ + 4)+0(^), (26) 



5ZZ = ~Cf^ — -3 7s + 31n^ + 4 +«, (27) 



±n !L 

A " \ L U V "^c / 

where fi is the renormalization scale, 7^ is the Euler's constant, and Cp = | for SU(3) C . 

Among the eight diagrams in Fig. (2j there are two self-energy diagrams, four triangle dia- 
grams and two box diagrams. The divergences of these diagrams can be analyzed as follows: 
the self-energy and the triangle diagrams contain UV divergences, while the box diagrams 
possess IR divergences. In addition, Coulomb singularities may arise from the box diagrams, 
due to the exchange of longitudinal gluon between c and c. Our calculations indicate that 
the UV divergences from the self-energy and the triangle diagrams are canceled by that from 
the renormalization constants of the quark field and the mass and IR divergences from the 
box diagrams are canceled by that in 5Zq. As one will see, the remaining Coulomb singu- 
larities will be canceled by Coulomb singularities in NRQCD matrix elements by matching 



condition. In our practical calculations, the 'Feycalc' package 12[ has been employed to take 



the trace of the 7— matrix and to reduce tensor integrals into scalar ones. Furthermore, we 



carry out the IR-safe scalar integrals by using G.'t Hooft and M. Veltman's method 13[. For 



the IR-divergent scalar integrals, we take the results from Ref. 



. 1 



After canceling all these 



UV and IR divergences, we present the NLO cross section of e + e —> cc( x 5*o) + 7: 

o = a (0) | q=0 (l + -C F ^ + *C(r)) , (28) 

\ V TV 71 J 

where a^lq^ is the LO cross section (fT9|) with q = and C(r) is a function of r and is 
both UV and IR finite. The explicit expression for C(r) is given in Appendix lAl 
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Subsequently, we turn to deal with the same process in the NRQCD factorization for- 
mulism. The main ingredient of this calculation is to obtain the NRQCD matrix element to 
the NLO in a s , which reads [lj], 

|(cc(%)|^ x |0)| 2 = 24m 2 (l + -C F ^) , (29) 

\ V TV J 

where the Coulomb singularity term emerges from the Coulomb gluon exchanging between 
the quark and the anti-quark. The total cross section of e + e~ — > cc^So) +7 in the NRQCD 
factorization formulism can be readily expressed as 

a = 24F l ( 1 S )(l + -CA). (30) 

\ V TV J 

By equaling (128jl with (1301 . we determine the NLO short-distance coefficient: 

F 1 ( 1 S ) = F}°\ 1 S )(l + ^C(rj), (31) 

from which we see that the Coulomb singularity cancels in the short-distance coefficient as 
expected. 1 + — C(r) is usually referred to the K— factor. 

Inserting this coefficient into ([2]), we reach the final result for the cross section to the 
NLO QCD correction in a s and the LO in v 2 : 



a = a (°) {1 + ^L C (r) + 0(a 2 s )), (32) 

7T 



where a*- ) is given by ( 123]) with q = 0. 

It is instructive to look at the behaviors of C(r) in two limits of r — > and r — > 00. 
The limit of r — > corresponds to a very high energy production of t] c . In this limit, C(r) 
reduces to 

lim C(r) = -- r(9-61og2)logr + 9(3 + log 2 2-31og2) + 7r 2 l . (33) 

r— >0 9 



The coefficient of logr term agrees with that obtained in Ref. [15| . 

The limit of r — 00 corresponds to the NLO QCD corrections to the process of r] c decay 
into two photons. To achieve this limit, one needs to recalculate the K— factor for r > 1. 
We define C'(r) in this region via 



° {r ^ VQl \ = 1 + ^C'(r) 



cr(°)(r* t)qj) ' n w v n ^ 
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whose expression is given in Appendix [A] With the expression of C'(r) in hand, we take the 
limit of r — > oo for C'(r), which reads 

7T 2 

lim C'(r) = CV(-5 + — ). (35) 

r— >oo 4 

It is encouraging that (I35I) is consistent with the QCD correction to the cross section of 
77- Comparing our results with that given in Ref. [4|, we found that the expressions 
given in their Eqs. (23) (25), from which they obtained the two limits of r — > and r — > oo, 
were incorrect. Actually, it is easy to check that the results given in their Eqs. (23) (25) go 
to infinity by taking these two limits. Therefore we claim that there might be some typos 
in their results. 



IV. e+e ->■ XcJ + lhard 

In this section, we apply the NRQCD factorization formulism ([2]) to the processes e + e~ — > 
Xcj + lhard with J = 0, 1, 2 and determine the short- distance coefficients Fi( 3 Pj) to the NLO 
in a s . As mentioned in Sec.HIl these coefficients can be determined by matching the processes 
e + + e~ — > cc( 3 Pj) + 7 with the cc near the threshold and in the color-singlet 3 Pj states. 



A. The LO short-distance coefficients in a s 



We first determine the LO coefficients in a s by calculating the cross sections of 
e+ G°ei) + e~(p e2 ) -» c(pi)c(p 2 )( 3 Pj) + j(k) both in the QCD and in the NRQCD factoriza- 
tion formulism. Two tree- level diagrams contributing to the processes e + (p e i) + e~(p e 2) — > 
c(pi)c(j)2)( 3 Pj) + 7(fc) are illustrated in Fig. [TJ To carry out the calculations, we need to 
construct the polarization vectors and tensors for \i an d X2 respectively. Three polarization 
vectors e^ Q ±1 j for xi can be explicitly expressed as(provided the hadron moving along z 
axis.) 

^ = -1(0,-1,-1,0), e^ = -l(0,+l,-z,0), e^ = ^-(|P|,0,0,E P ). (36) 
Then five polarization tensors for X2 are readily written in terms of e^ ±1 -, as 

t ±2 — fc ±l fc ±l' fc ±l — V fc ±l fc 1 fc fc ±l^ 

% = 7^(^-1 + 2^ + 6^). (37) 
12 



The amplitudes of e + e~ — > cc( 3 Pj) + 7 can be picked out from (fT3|) by extracting the 
contributions corresponding to the cc( 3 Pj) production. To this end, one needs to expand 
the amplitudes to the LO in q and to project it into the 3 Pj states with helicity \ H . It can 



be done by using the formulae 



id: 



A^{\ H ) = |q| l-Tr[A^U lv }V^ , (38) 
where II^ is defined via Hi = I\.\e u and the projection operators Vj U are defined by 

1 pn pv 

— —f—n^ + ———') 

Aivpa p * 

1 ~ 2V2E PXhC " 

V? = e m g. (39) 

With these formulae, it is straightforward to calculate the cross sections of the processes 
e + e~ — > cc( 3 Pj(Xh) + 7 with a specific helicity Xh for cc. The differential cross sections of 
these processes in the LO in a s yield 

da® A8e%a 3 n 2 



dx s 2 m 2 (l — r) 
where the functions Fj(\h) are 



W 2 Fj(\ H ), (40) 



Fo(0) = ^-^(1 + * 2 ) 
Fi(0) = (1 



.2> 



Fi(±l) = r(l -x%, 

3 



F 2 (0) = ^(1+x 2 ) 



,2> 



F 2 (±l) = r(l - x\. 
F 2 {±2) = r 2 (l + x 2 ). (41) 

Since relativistic corrections of 3 Pj production will not be calculated, we have set Mh = 
2E = 2m. Our formulae fj4*0l agree with that obtained in Ref. [rsl]. 

We now turn to evaluate the cross sections of the same processes in the NRQCD factor- 
ization formulism (T5|). The matrix elements for cc( 3 Pj) production are given by 

(cc( 3 Pj(\ h ))\^(}Cs Pj ■ ej) X \0) = 2V6m|q|. (42) 
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Thus, the differential cross sections over x within the framework of NRQCD factorization 
formulism read 

^ (0) -^(0),s 



dx 

,(0)/ 3 ; 



24Fr( 3 Pj(A„))|q| 2 , (43) 



where F{ ( Pj(Xh)) denote the LO short-distance coefficients. By equaling the cross sec- 
tions on the QCD side with that on the NRQCD side, we determine the short-distance 
coefficients: 

F?\*Pj(\ H )) = lif n \ Fj(X H ). (44) 



w?s 2 {l — r) 

Consequently, the differential cross sections of e + e~ — > Xcj(Xh)i in the LO in v and a s 
are written as: 

-fa- = 1 ^ 2 J) (0\x^Mxj(Xh)}(xj(Xh)\^1Cjx\0}, (45) 

where the xj states in the matrix element are normalized relativistically. To express our for- 
mulae as the product of short-distance coefficients and standard non-relativistic normalized 
matrix elements, we can use the relations between relativistic normalized matrix elements 
and non-relativistic normalized matrix elements jlo| 

(0|xt/Cs Pj V|if) = V / 2Mfl : (0|x t /C3p J V|iy) BBi) (46) 

where \H)bbl is the non-relativistic normalized matrix element used in Ref. 1] and M# is 
the mass of H. In the LO in v, we set Mh = 2E\ q =Q = 2m. 

The differential cross sections of e + e~ — > Xcj(Xh) + 7 can be reexpressed as: 

da® 8e A a 3 n 2 ~ 

= m3^(l-r) Fj(Ag)(0|X ^ IXcj)(Xcj(AH)l ^ JX(Ag)|0) ^ (47) 

The total cross sections are achieved by integrating over x and summing over all the 
helicity states Xh of the xj- We present the final results: 

{0\x^ji J \Xcj){Xcj\i'^jX\^)BBL , (48) 

where Fj for J = 0, 1, 2 are given by 

(0) 32e 4 Q a 3 7r 2 (l-3r) 2 

" 9m¥(l-r) ' 
f(0) = Ue%a^{l + r) 

1 3m 3 s 2 (l — r) ' ' 

m 64eSa 3 vr 2 (l + 3r + 6r 2 ) 
F 2 (0) = 5 -. 49 

2 9m 3 s 2 1 -r) 



14 



B. Short-distance coefficients in the NLO in a s 



In this subsection, we determine the short- distance coefficients in the NRQCD factoriza- 
tion formulism for processes e + e~ — > Xcj + 7 up to the NLO in a s . As done in the case 
of r\ c production, we calculate the cross sections of the processes e + e~ — > cc( 3 Pj) + 7 both 
in the QCD and in the NRQCD factorization formulism to one-loop level. Feynman dia- 
grams responsible for the processes are the same as that of 1 So production as illustrated in 
Fig. |2j We still carry the renormalization in on- mass-shell (OS) scheme. The wave function 
renormalization constant and mass renormalization constant have been given in ( 126]) f )27|) . 

The UV and IR divergences of each diagram are similar with that of rj c production. The 
self-energy and triangle diagrams contain ultraviolet (UV) divergences. The box diagrams 
have both infrared(IR) divergences and Coulomb singularities. Combining contributions 
from the eight diagrams and the counterterms, all divergences cancel except Coulomb sin- 
gularities which will be canceled by Coulomb singularities in NRQCD matrix elements. The 
cross sections of the processes e + e~ — > cc( 3 Pj, Xh) + 7 in the NLO in a s are given by 

da 48e> 3 vr 2 ~ / a a \ 

|q| Fj(\ H ) H C F h — Cj H (r) }, (50) 



dx m 2 s 2 (l — r) ' \ v it ir 
where the index A# denotes the helicity of Xcj, Fj(Xh) have been given in previous subsection 
and the expressions for Cj H (r) are listed in Appendix El 

One can immediately get the total cross sections by integrating over x and summing over 
the polarizations of Xcj- As a result, the cross sections yield 



7r 2 „ a, a, 



a = FT |q| 2 1 + ~Cf- + -Cj(r) (51) 

V IT 71 



where Cj(r) for J = 0, 1, 2 are given by 

C (r) = C° (r), 

C?(r)+rCKr) 

C°(r) + 3rCl(r) + 6r 2 Ci(r) 
C2[r} ~ l + 3r + 6r 2 " (52) 

We now calculate the cross sections of the same processes in the NRQCD factorization 

formulism. The NRQCD matrix elements of P-wave operators in the LO in v and NLO in 

a s are given by 

|(cc( 3 Po)|/C3 Pj |0)| 2 = 24m 2 |q| 2 (l + —C F ^) . (53) 



V TT 
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With it, one readily obtains the cross sections of e + e — > cc( 3 Pj) + 7 on the NRQCD side: 

a = 24F 1 ( 3 P J )|q| 2 [l + -CV^]. (54) 



7T 



The short-distance coefficients Fi( 3 Pj) can be determined by equaling (l5Tj) with (154|) . 

F^Pj) = Fi°\ 3 Pj)[l + ^Cj(r)]. (55) 

7T 

Finally, the cross sections of the processes e + e~ — > Xcj + Ihard in the NLO in a s can be 
expressed as 

°{Xcj) = a {0) ( Xc j)(l + ^Cj(r)), (56) 

where <rf] denote the LO cross sections given in fj4*8|) . 

It is interesting to look at the behaviors of Cj(r) in two limits of r — > and r — > 00. The 
limit of r — » is related to the XcJ production in a very high energy. In this limit, we find 

2 

limCo(r) = -- [(3 - 61og2) logr + 91og 2 2 - 331og2 + tt 2 1 , 

r-s-0 9 

limCi(r) = -\ [(9 -6 log 2) logr + 9 log 2 2- 15 log2 + vr 2 + 2ll , 

r->0 9 

hm C 2 (r) = -- [(3-61og2)logr + 91og 2 2 + 31og2 + 7r 2 + 18l . (57) 

r— s-0 9 

The limit of r — > 00 is related to the NLO QCD corrections to the cross sections of Xcj 77- 
As in the case of r] c , one needs to recalculate the K— factor for r > 1. We define C' cJ {r) via 

4^4 = l + ^(r). (58) 

Explicit expressions of C' cJ (r) are listed in Appendix [Bl In the limit of r — > 00, Cj(r) go to: 

r— too 3 4 

lim C' 2 = -AC F . (59) 



r— >-oo 



The formulae (1591) are desirable, since these are nothing but the QCD corrections to the 
processes Xo,2 - > 77- In (|59l - we did not take the limit of r — > 00 for C{(r), since the cross 
section of Xci decay into 77 vanishes due to the Yang's theory. 
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V. e + e ^ rib(xbj) + Jsoft 

As discussed in Sec. II, for bottomonium production at the B factories through process 
e + + e~ — » bottomonium +7, the emitted photon is soft comparing to the b quark mass. The 
cross section of the process is expressed as in the framework of the NRQCD factorization 
formulism. The short- distance coefficients -F^ 3 ^) and G2{ 3 S\) in © can be determined by 
matching the cross section of the process e + (p e i) + e~{p e 2) — > c(pi)c(p2) with the cc in the 



color-singlet 3 Si state predicted by QCD and by the NRQCD factorization 
has been calculated up to order a 2 . Here we list only the -F 2 ( 3 S'i) result j^jj^] 



ormu 



|7j Q|20j: 



ism. It 



An 3 e 2 n a 2 5(s - Am 2 ) ( a<1 IAOtt 2 2ml a 2 

F 2 3 5i = Q ^ L)i_ ACf ^+ -117.46 + 0.82^ + ^— log — PI 

m I 7r 27 n J n z 

(60) 

The matrix elements in (|9]) describe probabilities of a point- like 3 Si bb state transiting into 
the 1 S'o or the 3 Pj state through a magnetic Ml transition or an electric El transition 
by emitting one soft photon. Although they are QED processes, the matrix elements are 
non-perturbative in nature. One may evaluate those Ml and El transition processes by 
perturbative QCD and relate the matrix elements appearing in (Q to those appearing in 
(j2J). However, this kind of calculations is not reliable. Some additional Coulomb singularities 
may arise from gluon exchange between almost on-shell b and b quarks. Terms like a s /v 
appearing in radiative corrections spoil fixed order perturbation calculations on those matrix 
elements. Even at tree level, it is proportional to a factor s/4m 2 — 1 for the 77^ production and 
l/(s/4m 2 — 1) for the rj c j production. Thus the results are very sensitive to the quark mass. 
To gain a reliable prediction on the matrix in (Q, one need to invoke some nonperturbative 
methods such as lattice QCD. 



VI. NUMERICAL RESULTS AND DISCUSSION 



In this section, we carry our numerical calculations on the production cross sections of 
e + e~ — > H +7 (H = r] c , rj c (2S), Xcj) at the B factories using the results derived in section III 
and section IV. To this end, we need to determine some relevant parameters such as quark 
mass, ot s (/j,), and those NRQCD matrix elements. We will estimate them below. 
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A. Input parameters 



At the B factories, we take \fs = 10.6GeV and a(fi = y/s) = 1/131. The quark mass is 
taken as m c = 1.4GeV. There is an ambiguity in choosing the value of [i in a s since there are 
several energy scales involved, such as s/s and heavy quark mass. In the limit of r — > 0, one 
can use the leading log approximation (LLA) to resume those large logarithmic terms such as 
(a s (-\/i) In v^i/m)™^, 15 1. This accounts for the effects of the running of the strong coupling 
constant a s from energy scale to m. However, in doing this approximation, one has 
omitted constant terms and terms suppressed by power of r. To recover these contributions, 
one needs to carry out calculations to higher loops and to include contributions from higher 
dimensional operators. Those are much more complicated. In our full one-loop results, we 
keep all these terms but leave an ambiguity in choosing the value of /i in a s (fi). To estimate 
the uncertainty caused by it, we take three different values of /i. They are a/s, 2m c , and m c . 
The corresponding values of a s (fi) are a s (y/s) = 0.17, a s {2m c ) = 0.24, and a^m^] = 0.30. 



As for the hadron matrix elements, we will take them as same as those used in 



(o\x^\vc)(vc\^x\o) 


= 0.437GeV 3 , 


(0\xty\Tk&S))(ric(2S)\tfx\0) 


= 0.274GeV 3 , 


i|( Xc0 |^(-lD.a)x|0)| 2 


= 0.051GeV 5 


l|(x cl |^(_lB X( x. e// ) X |0>| 2 


= 0.060GeV 5 , 


E^|^(-^V)e|) X |0)| 2 


= 0.068GeV 5 . 



(61) 



For the r) c production, in order to assess the contribution from the relativistic correction, 
we will evaluate the NRQCD matrix element corresponding to the relativistic correction by 



using the Gremm-Kapustin 



23[ relation, 



Re(0|x t (-|D)V|7 7c )(r ?c |^ t x|0) = mV Re(0\ X ^\Vc)(Vc\^x\0)- 



(62) 



We take v 2 = 0.13 for rj c as that in Ref. 24j and the correspondent NLO NRQCD matrix 
element yields: 



Re(0|x T (-2 D )^l^)(^|^x|0) 



O.lllGeV . 



(63) 



For the r] c (2S) production, the LO NRQCD matrix element, which is decided by fitting 
data from T[ijj(2S) —> e + e~] based on heavy-quark spin symmetry, is accurate only in the 
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TABLE I: Predicted cross sections (in fb) for the production of r/ c , rf c (2S), and XcJ with various \x. 
cr(°) and a are the cross sections with the LO and the NLO short-distance coefficients, repectively. 
-joy — 1 is relative corrections caused by the NLO corrections. 



Li = yfs /i = 2m li = m 





_2 1 


a 


_2 1 


a 




a 


Vc 


83.3 


-13.0% 


72.5 


-18.4% 


68.0 


-23.0% 


64.1 


Vc(2S) 


52.2 


-13.0% 


45.4 


-18.4% 


42.6 


-23.0% 


40.2 


XcO 


1.19 


10.4% 


1.31 


14.7% 


1.36 


18.4% 


1.41 


Xcl 


14.3 


-17.1% 


11.9 


-24.1% 


10.9 


-30.1% 


10.0 


Xc2 


6.28 


-48.9% 


3.21 


-69.0% 


1.95 


-86.2% 


0.87 



LO in v as pointed out in Ref. [2|. Therefore we will not consider the relativistic correction 
to T] C (2S) production for consistence. 

In the following subsection, we will use these input parameters to calculate the production 
cross sections. 

B. Numerical calculations 

With the parameters given above, we can evaluate the cross sections for the production 
of rj c , r] c (2S), and Xcj by using Eqs. (123]) . (|32|) . ( 148|) . and (!56l) . For comparison, we list the 
results of the LO and the NLO in a s in Table [H For the NLO results, we have taken \x in 
a s (n) to be y/s, 2m c , and m c . In TabelHJ we have omitted relativistic corrections. From 
the table, we see that most of the QCD radiative corrections are negative and considerable. 
For some processes like \2 production, the corrections are unexpectedly large. 

Then we can estimate the relativistic corrections by Eq. (1231) . After accounting for the 
NLO radiative QCD correction(we take [i = 2m c in a s (/i)) and the relativistic correction, the 
total cross section of the rj c production decrease from 83.3fb to 58.2fb, where the NLO radia- 
tive QCD correction and relativistic correction contribute -15.3fb and -9.81fb respectively. 
We see that both of them give considerably negative contributions. 

These theoretical predictions for the cross sections receive uncertainties from the am- 
biguities of the values of input parameters. These parameters are the c quark pole mass, 
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TABLE II: Predicted cross sections (in fb) for the production of rj c , r/ c (2S), XcJ by including 
theoretical uncertainties. 



Vc 


77,(25) 


XcO 


Xcl 


Xc2 


68.0^-2 




1 Sfi +0 - 26 




1 Qc+1-85 
± - yo -1.56 



the NRQCD matrix elements, and the \x in ct s (/j,). The c quark pole mass may effect both 
the short-distance coefficients and NRQCD matrix elements. Especially, the matrix element 
of the relativistic correction term is very sensitive to the value of the pole mass when the 
Gremm-Kapustin|23 relation ( |62|) is used to estimate its value. For more details, we refer 
the reader to Ref. [2j. Another theoretical uncertainties arise from the higher order QCD 
corrections and higher order relativistic corrections. Including these corrections, the fi de- 
pendence and uncertainties of the matrix elements may be reduced. Varying the values of 
all these parameters in reasonable ranges and omitting higher order correction terms both 
in a s and in v, we can estimate theoretical uncertainties on the predicted cross sections. 
The results are listed in Table [TTJ 

In e + e~ collisions, the heavy quarkonium with C— parity even can also be produced via 
the double photon processes. At very high energy limit, this process can be described by the 
equivalent photon approximation (EPA). Comparing to the processes e + e~ — > H + 7 (H = 
r] c ,r] c (2S),Xcj), the total production cross sections of the double photon processes are sup- 
pressed by additional power of a but enhanced by log 2 s/m 2 e [25| . The cross sections of 
the s— channel processes e + e" — >■ H + 7 (H = r] c ,r] c (2S),Xcj) are suppressed by M\j s 2 . 
Combining all these factors, at the B factories, the production cross sections of the heavy 
quarkonium with C— parity even via the double photon process are larger than that via asso- 
ciated production with photon. However, the signals in detectors from these two production 
mechanisms are very different. Especially, the produced heavy quarkonia via double photon 
processes are predominated by small Pt events at high energy which are hard to be detected 
from their decay products. On the other hand, the produced heavy quarkonia via associated 
production with photon are insensitive to the Pt- Moreover, the events can be reconstructed 
by measuring the photon assuming the photon detector is good enough. A thorough analysis 
on these processes in the framework of the NRQCD factorization formalism will be presented 
elsewhere. 
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Appendix A: The NLO radiative corrections to short distance coefficients for the 
rj c production 

The definition for C(r) appeared in Eq. fl28l) : 

2[30r 2 -(84 + 7r 2 )r + 27r 2 + 54] 8(2r - 3) , .2 . 4, .1-6. 
CW = " J 9(2-r)(l-r) ' + W^f ~ 2) " b ^ 

+ 3(^1) [(1 + bg2( TT*> " bg2( ' - 1)] + W^) U ^), (Al) 

where we use r = b = a/1 — 4m 2 /s and m = m c in Eq. ( 1B1I) . In addition, 

dilogarithm(Li 2 (x)) is used in the Appendices lAl and IB1 which is defined as 

Li 2 = - / -log(l-z), (A2) 
Jo z 

where x < 1. 

The definition for C"(r) appeared in Eq. (134"]) : 

7r 2 (r-2)(3r + 4)-12(r-l)(5r-9) 2(2r - 3) 2 

C (r) " ° F 12(r-l)(r-2) + ° F J^W g( " r } 

CV r _ ,2 



{Li 2 (- - 1) - tan" 1 (_L=)[6^~1) 



r - 1 r v /r - 1 

+ (r + 2)tan- 1 (^i=)]}, (A3) 
— 1 
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Appendix B: The NLO radiative corrections to short distance coefficients for the 
Xj production 



1- Xo 



C° (r) 



7r 2 (4r 4 - 7r 3 + llr 2 - r + 1) 



+ 



+ 



9(3r 3 -r 2 + 3r-l) [ (1 - 

6r (7r 4 - 25r 3 + 31r 2 - 25r + 24) 
(2-r) 2 

1 f2(4r 4 -7r 3 + llr 2 -r + l)log 2 (f -1) 



3(r-l) 2 (3r-l) 



1 + r 2 



1 — 6 4(4 — r") 
+ (3r 3 - 6r 2 - 3r + 2) log 2 (p-^) + — ^(3r 2 - 12r + 5) log( 

8(2r 4 - 5r 3 + 6r 2 - 2r + 1) log(f - 2) log(f - 1) 
(1 +r 2 ) 

Br 5 - 50r 4 + 125r 3 - 152r 2 + 87r - 18) log(^ - 2) 



1-^ 
1 + b' 



+ 



+ 



(r - 2) 3 

4 [(4r 4 - 7r 3 + llr 2 - r + 1) Li 2 (^) - 2 (2r 4 - 5r 3 + 6r 2 - 2r + 1) Li 2 (2 - 

3(3r - 1)(1 -r) 2 (l + r 2 ) ' 



2. Xi 



Cf\r) 



+ 
+ 
+ 



7T 2 (-4r 2 + 19r - 5) - 36(r - l) 2 (-3r 2 + 3r - 2) log 2 (|g) 

9(r-l) 2 + 3(r-l) 2 

4[(r - 2)(2r 2 - r + 1) log + b (2r 2 - 5r + 1) log (2 - 2)] 

36(r - 2)(r - 1) 

[(4r 2 - 19r + 5) log (f - l) - 2 (r 2 - 16r + 3) log (f - 2)] log (f - l) 

3(r- l) 2 

2 (4r 2 - 19r + 5) Li 2 (^) - 2 (r 2 - 16r + 3) Li 2 (2 - fj 

3(r- l) 2 ' 
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2[n 2 (r - 2)(r 2 - r + 2) - 6(r - l)(4r - 7)] 4(r 2 - 2r - 1) log (±=§) 



9(r - l)(r - 2) 



36(1 -r) 



+ 



8(3r 3 - 10r 2 + 9r - 1) logg - 2) (r - 2)(r 2 + 1) log 2 (±^) 

3(r-2) 2 (r- 1) + 3(r - l) 2 

2 log (| - 1) [2(1 - 2r) log(f - 2) + (r - l)(r 2 - r + 2) log(f - 1)] 

3(r- l) 2 

4[(1 - 2r)Li 2 (2 - \ ) + (r - l)(r 2 - r + 2)Li 2 (^)] 
3(r-l) 2 



(B2) 



3- X2 



2r 2 [24(r-l) 2 + 7r 2 (r + 3)] 8r 2 [4Mog (| - 2) - (r - 5) log (|g)] 
9(r-l) 2 + 36(r-l) 

2r 2 {log (2 - 1) [(r + 3) log g - l) - 2(r - 1) log (3-2)] - 41og 2 (jg )} 

3(r- l) 2 

4r 2 [(r + 3)Li 2 fe) - (r - l)Li 2 (2 - \ )] 
3(r-l) 2 



tt 2 (2 - r)(2r 2 - 20r + 3) - 12r(r - l) 2 (3r 2 + 21r - 2) log 2 (g) 





9(r - 2)(r - 


-l) 2 3(r-l) 2 




log (? - 


1) [(2r 2 - 20r 


+ 3) log (2-1) -2(5r 2 + r + l)log(f 


-2)] 






3(r- l) 2 




4[(19r 4 


-3)(r-2) 2 log( 


;i=|) + 6(10r 3 - 49r 2 + 51r + 10) log(f ■ 


-2)] 






3b(r - 2) 2 (r - 1) 




2[(5r 2 + 


r + l)Li 2 (2- 


f ) - (2r 2 - 20r + 3)Li 2 fc)] 





3(r-l) 2 
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7T 2 (2r 4 + 9r 3 + 27r 2 - llr + 5) (r - 2) 2 + 12(r - l) 2 (2r 3 - llr 2 + 6r + 12) 



C 9 u (r 



2V ; Q/V _ _ 1 \2 



9(r-2) 2 (r- l) 2 

(3r 3 + 30r 2 + 3r + 2) log 2 (±=g) (2r 4 + 9r 3 + 27r 2 - llr + 5) log 2 (2 - 1) 

3(r- l) 2 + 3(r- l) 2 

8(6r 5 - 38r 4 + 60r 3 + 16r 2 - 69r + 6) log (f - 2) 4(21r 2 + 18r - 1) log(^|) 



3(r-2) 3 (r-l) 36(r - 1) 

2 (2r 4 + 9r 3 - 9r 2 - llr + 3) log (f - 2) log (f - l) 
3(r- l) 2 

2 2 
[(-2r 4 - 9r 3 + 9r 2 + llr - 3)Li 2 (2 - -) 



3(r - 1) 



+ (2r 4 + 9r 3 + 27r 2 - llr + 5)Li 2 (- — -)]. (B3) 



4. the expressions for C'Ar) 



For C' (r): 

C'Jr) = — — ° F -[-6r(7r 2 -25r + 24)(r- l) 2 

0V ; 6(r-2) 2 (r-l) 2 (3r- 1) 1 v A 1 

+ vr 2 (r - 2) 2 (3r - 1) - 6vr(r - 2)Vr - l(3r 2 - 12r + 5)] 

2C F (2r - 3)(4r 3 - 15r 2 + 19r - 6) log(2 - 2) C F U 2 (f - l) 



+ 



(r - 2) 3 (r - l)(3r - 1) (r - If 



C F 



+ (r _ 1)2( 3 r _ 1} tan- 1 ( v / ^T)[2v / ^T(3r 2 - 12r + 5) 



+ vr(3r 3 - 6r 2 - 3r + 2) - (3r 3 - 6r 2 - 3r + 2) tan" 1 (y/r - l)]. 



For C[{r): 



C'Ar) = — x , Cf — r[-12(r- l) 2 (3r 2 - 2r - 7) 

1V ; 12(r-2)(r-l) 2 (r + l) L V ; V ; 

+ vr 2 (r - 2)(r 3 + r 2 + 2) + 12vy/r - l(r - 2)(2r 3 - 2r 2 + 3r + 1)] 

C F (2r 4 - 3r 3 - 9r 2 + 16r - 2) log (2 - 2) C F (r 3 + r 2 + 2) Li 2 (2 - l) 



(B4) 



( r _ 2 )2( r 2_i) 2(r-l) 2 (r + l) 



C F 



+ 7 7^7 7T t a n- 1 (v / ^I)[-2v / ^ r T(2r 3 -2r 2 + 3r + l) 

(r — l) z {r + 1) 



7T 



(2r 3 - r 2 + r + 2) + (2r 3 - r 2 + r + 2) tan" 1 ( v / r - 1)] . (B5) 
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For C' 2 {r): 
C' 2 (r) = 



-c t 



— [12 (24r 4 - 91r 3 + 79r 2 + 6r + 12) (r - 1] 



12(r - 2) 2 (r - l) 2 (6r 2 + 3r 
+ 7r 2 (r - 2) 2 (9r 3 - 21r 2 - 6r - 2) + 12?r\/r - l(r - 2) 2 (12r 3 - 24r 2 - 9r + 1)] 
C F (30r 5 - 157r 4 + 249r 3 - 76r 2 - 78r + 12) log (2 - f ) 



(r -2) 3 (r- l)(6r 2 + 3r + l) 
C F (9r 3 - 21r 2 - 6r - 2) Li 2 (f - l) 
2(r - l) 2 (6r 2 + 3r + 1) 
C 

— tan _1 (Vr - l)[-2>/r - l(12r 3 - 24r 2 - 9r + 1) 



(r- l) 2 (6r 2 + 3r + l) 
+ ?r(-6r 3 + 15r 2 + 9r + 2) + (6r 3 - 15r 2 - 9r - 2) tan -1 (Vr - 1)] 



(B6) 
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